Notes On The Post-Newtonian Limit Of Massive Brans-Dicke Theory 
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We consider the Post-Newtonian limit of massive Brans-Dicke theory and we make some notes 
about the Post-Newtonian limit of the case ui = 0. This case is dynamically equivalent to the 
metric f(R) theory. It is known that this theory can be compatible with the solar system tests if 
Chameleon mechanism occurs. Also, it is known that this mechanism is because of the non-linearity 
in the field equations produced by the largeness of the local curvature relative to the background 
curvature. Thus, the linearization of the field equations breaks down. On the other hand we know 
that Chameleon mechanism exists when a coupling between the matter and the scalar field exists. 
In the Jordan frame of Brans-Dicke theory, we have not such a coupling. But in the Einstein frame 
this theory behaves like a Chameleon scalar field. By confining ourselves to the case u: = 0, we 
show that "Chameleon-like" behavior can exist also in the Jordan frame but it has an important 
difference compared with the Chameleon mechanism. Also we show that the conditions which lead 
to the existence of "Chameleon-like" mechanism are consistent with the conditions in the Post- 
Newtonian limit which correspond to a heavy scalar filed at the cosmological scale and a small 
effective cosmological constant. Thus, one can linearize field equations to the Post-Newtonian order 
and this linearization has not any contradiction with the existence of " Chameleon- like" behavior. 

PACS numbers: 98.80.Es,98.65.Dx,98.62.Sb 
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I. INTRODUCTION 

The dynamics of massive Brans-Dicke (BD) theory can 
be determined by the following action 



S= 2k^J dAx ^~9^R - -jd^f - V (</))] + S m [g, V] 

(1) 

where k 2 = 8irG and S m is the matter action. In the 
original BD theory, the mass term V (4>) is zero and so 
the scalar field is massless. The Post-Newtonian (PN) 
limit of massless BD theory has been investigated in 
It is known that the massless BD theory can not always 
produce the corresponding GR case in the limit lu — > oo. 
And there exist some exact solutions that doesn't go over 
to the Einstein's theory in this limit Banerjee et al Q 
argued that this is due to the zero trace of the energy- 
momentum tensor of matter. But even for a non zero 
trace of energy-momentum tensor, the limit of uj — > oo 
doesn't lead to GR necessarily Q ■ Whether the BD the- 
ory behaves like GR at oj — > oo or not, depends on valid- 
ity of the main assumption of PN formalism. For exam- 
ple, in the PN limit we assume that the perturbation of 
the scalar field due to the local gravitating system under 
consideration is very small compared to its cosmological 
background value. If this assumption fails then this the- 
ory will not behave as GR in the PN limit. The situation 
is more complex in the massive BD theory even if the 
coupling constant uj is zero. The weak field limit of this 
case is the subject of controversy. More specifically, sev- 
eral authors claimed that the metric f(R) theory (which 
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is equivalent to the massive BD theory with u) = 0) is 
characterized by an ill-defined Newtonian regime Q . 

It is known that, the Jordan frame of BD theory is 
related to the Einstein frame via a conformal transfor- 
mation. In the Einstein frame there exists a conformal 
coupling between the matter and the scalar field which 
leads to an interaction between them 0. Hence, this 
theory in the Einstein frame can be considered as an in- 
teracting quintessence modelQ . On the other hand, this 
interaction leads to a density-dependent scalar field mass 
i.e. the associated mass of the scalar field can change 
with environment Q. This behavior is the reason for 
naming such a theory the Chameleon theory Q. The 
mass of the scalar field can be very large at the dense 
places and so the scalar field interaction with the matter 
can be strongly suppressed. Also, under the condition 
named the thin-shell condition Q the scalar field out- 
side the source is produced only with a very thin shell of 
matter near the surface of the source. In the other words, 
the interior part of the source has no contribution to the 
generating of the scalar field. 

An important question can be arisen here. We know 
that some ambiguity exist for the physical equivalence of 
the Einstein frame and the Jordan frame. Can we con- 
clude the existence of the Chameleon mechanism in the 
Jordan frame from its existence in the Einstein frame? 
We expect that one should see this behavior directly 
in the Jordan frame without switching to the Einstein 
frame. Also it is natural to expect the trace of this be- 
havior in the PN limit. However, it is claimed in the 
literature that when Chameleon mechanism occurs, lin- 
earization of the field equations (which is necessary in the 
PN limit) breaks down and the behavior of the scalar field 
is governed by non-linear dynamics @. But, by working 
directly in the Jordan frame we show that there is no 
Chameleon mechanism in the form of 0] but a differ- 
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ent version of it which wc named it "Chameleon-like". 
We verify that linearization of the field equations to the 
PN order has not any contradiction with the existence of 
the "Chameleon-like" mechanism and the required con- 
ditions for the existence of this mechanism are consistent 
with the required conditions in the PN limit for the via- 
bility in the solar system. 

The outline of this paper is as follows: In Sec. II we 
review the PN limit (expansion) of the massive BD the- 
ory and we look for non-linear terms responsible for the 
Chameleon mechanism. Although, we will show that 
such terms do not exist. In Sec. Ill wc confine ourselves 
to oj = case and show that " Chameleon-like" behavior 
can exist in this theory. Finally, we compare the result 
of PN limit and Chameleon-like behavior and we show 
that they are consistent with each other. Also, we clarify 
the range of validity of the PN expansion and show that 
this expansion dose not break down when the Chameleon 
mechanism occurs. 
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The full derivation of the PN limit of the massive BD 
theory can be found in [Toj]. Let us just review in brief 
the results of which are necessary for our work. The 
field equation (j4]) in the PN approximation is 



[V 2 -m 2 ] <pW{x t t) 



k 2 p 
' 3 + 2uj ' 



(6) 



Where 



is the mass associated with the scalar field at 



the cosmological scales given by 

2 W ~ 



V 



2oj 



(7) 



in which prime denotes differentiation with respect to 



("2") 

In the Solar system the expression of h 00 and h\*' take 
the following forms far from the source 



,(2) 



II. POST-NEWTONIAN LIMIT OF MASSIVE 
BD THEORY 



In order to find some solutions of the field equations 
in the PN approximation, we expand the equations of 
motion around the background values of the metric and 
the scalar field. More specifically, we will use: = 
V»» + V. fl"" = *T -W V A (x, t) = (t) + ip (x, t) and 
V 0) =Vo + <pV ' + (f 2 V "/2 + .... It is worth to mention 
that we are working in the Post-Newtonian coordinate 
system Q. The PN limit of any scalar-tensor theory 
requires a knowledge of: goo to 0(4), goj to 0(3), % to 
0(2) and ip to 0(4). Note that in the PN limit we take 
v 2 ^ JJ ^ 0(2) where v is the characteristic velocity of 
particles and U is the Newtonian gravitational potential. 
Thus, we look for the solutions of the field equations in 
the form of a Taylor expansion as 



3oo 
9oj 
9ij - 



i + c 

h {3) 



h 



(4) 

00 ! 



k) + V 



(2) , 



(2) 



(4) 



Variation of the action (fTJ with respect to g^ v and 4> will 
yield to the following field equations respectively: 



. (2) _ r, GeffM Q Ae// 2 
"00 — z ~ 1 7, ' i 
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where A e /y = Vb/2</>o and 
k 2 
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Note that in the derivation of equation © we have used 
the assumption m 2 , > 0. Using the equation ©, the PPN 
parameter 7 is given by 



(2) 
Hi 
,(2) 



3 + 2lu - e" 
3 + 2u> + e~ 



(10) 



See also [111 ] in which a special potential has been investi- 



gated in the PN limit. If the scalar filed is very light then 
the parameter 7 is space- independent and as we expect, 
the massive BD theory behaves like massless BD theory 
i.e. this parameter takes the form 7 = (1 + oS)/{2 + lu). 

We are interested in the case u = which, as men- 
tioned before, is dynamically equivalent to the metric 
f(R) theory. In this case if the scalar field is very light 
then 7 = 0.5 which is in patent disagreement with obser- 
vation since 7 b ~ 1 [T2|. However, if 
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mlh 2 > 1 
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where L represents a typical experimental length scale 
(below the planetary scale), then the scalar field is heavy 
and its interaction is short-range and it will be hidden 
from the local experiments. In this case 7 is near to 
1. On the other hand, the cosmological constant term 
(Vo/6t/>o) r2 must be very small. Otherwise, this term 
will change the gravitational dynamics of local systems 
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such as solar system. Thus, for having an observation- 
ally acceptable theory (at least in the solar system) the 
following condition should be satisfied too 



V0 T2 



Li < l 



(12) 



where Li is a length scale of the same order or larger 
than the solar system. Thus, the massive BD theory with 
uj = can be compatible with the solar system tests if 
both conditions (fTT|) and (fT2"j) are satisfied. Note that 
this conclusion are also true for metric f(R) theory fl"o| . 

Some notes are in order here. It is claimed in the lit- 
erature that the Chameleon behavior is due to the non- 
linear terms in the field equations^. So, we have to 
reconsider the PN limit with more care. May be there 
exists some non-linear terms in the PN expansion which 
have been forgotten in the above considerations. The 
main assumptions on the scalar field in obtaining the PN 



limit are 



< 1 and 



.(2) 



-,(4) 



Let us 



<Po + <P K ' + V 
examine the first assumption. Can we use this assump- 
tion in the gravitational systems such as solar system? In 
fact, (p^> is the contribution of the scalar field due to the 
loc al g ravitating system. A comparison with the results 
of [13| may be useful here. In order to find the weak field 
limit of the metric f(R) theory, Chiba, T. et al [13| have 
expanded the Ricci scalar as R = Rq + Ri{r), where Rq 
is the background curvature and Ri(r) is the perturba- 
tion produced by a given source. In order to linearize the 
field equations they have used R\{r)/R <C 1. However, 
as we will discuss in the next section, this assumption 
is experimentally unacceptable near the Earth or at the 
solar system and therefore leads to a wrong weak field 
limit. Similarly, in the PN limit we linearize the field 
equation of (j> using the condition 
arbitrary system, p^ 



< 1. If, for an 



is the same order or greater than 
4>o then it is straightforward to show that the PN expan- 
sion fails and we will not able to use this approximation 
for that system. In this paper we do not want to find un- 
der which conditions this assumption is applicable, but 
we will show that when the thin-shell condition satisfied 
then 



<C 1 is also satisfied. Thus we expect that the 

physical results of Chameleon-like behavior and the PN 
limit be similar. 

Now, consider the second assumption. We know that 
the odd-order terms 0(1),0(3) can not exist in the g 00 
and gij component of metric. In fact, conservation of 
rest mass prevents terms of order O(l) and conservation 
of energy in the Newtonian limit prevents terms of order 
0(3). Appearance of other odd-order terms is theory de- 
pendent, for example, 0(5) can not exist in GR (because 
of conservation of energy in the PN limit) but 0(7) can 
exist @. Also, we find the appropriate expansion of the 
scalar field from its field equation. For example, consider 
the field equation of (j> field in the original BD theory 



U(j) 



8tt 



2w 
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(13) 



By expanding the RHS for a perfect fluid to the fourth 
order and using an appropriate gauge condition, we can 
easily show that only even-order terms can appear in the 
expansion of <p (i-e. <ft — ipo + f^+f^) [HI- However, in 
the massive BD theory the situation is not so trivial. In 
this case, consider equation (j4]) and assume the following 
expansion for the scalar field <j) 



^^„ + + ^ + ^ (3) + p (i} + 



(14) 



It is an easy job to verify that if n is even(odd) then Oip^ 
contains only even(odd)-order terms. Also, we know 



7<r 



h v h v 



+ 0(5) (15) 



and will use the following gauge conditions 

1 



a 0,v 



hi. 
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-K , 

2 V - J 



(16) 



Note that for our purpose it is not necessary to fix the 
explicit form of £Ro and However, we should keep 
in our mind that 9\ is of order 0(3) and IK, is of order 
0(2). Using these gauge conditions and arranging the 
terms with the same orders in we obtain 
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(17) 



Thus, it is obvious that the equation of motion of <j), un- 
like the masslcss BD theory, allows the odd-order terms. 
Also, it is interesting that there exists a nonlinear term 
in the equation of ip( 2 \ If ip^ ^ then this term will 
change the solutions of the metric components and also 
the scalar filed itself. But, we should take into account 
the metric field equations too. If one writes the 0-0 com- 
ponent of to the second order, then there will be a 
term with the first order which should be zero separately 



(V 2 + ^ > (1) = 



(18) 



Similarly, from i-j component of ([3]) we have didjtp^ 1 ' = 
and from its 0-j component dodjtp^ 1 ' =0. Thus, the met- 
ric field equations force ip^ to be zero. This conclusion is 
also true for tp^ and one can check it by writing equation 
© to the higher orders. As a final result of this section, 
one can be sure that the PN expansion obtained in [l(| 
is complete and there is not any non-linear term in the 
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field equations. So, we should see the trace of Chameleon 
mechanism in this PN limit. In the next section we show 
that "Chameleon-like" behavior can exist in this theory 
and then we show that how this mechanism and PN limit 
arc connected. 



III. CHAMELEON-LIKE BEHAVIOR IN THE 
MASSIVE BD THEORY WITH cj = 

In this section we show that the Chameleon-like be- 
havior can exist in the massive BD theory with vanishing 
coupling constant. Our main purpose is to address the 
question that: Can we see the Chameleon mechanism in 
the PN limit? 

As we discussed in the introduction, for Chameleon 
scalar fields the associated mass is dependent to the mat- 
ter density of the environment and so the scalar field is 
short-range in the dense places and is long-range in the 
low densities such as cosmos Q. The main reason for 
such behavior is that in these theories the scalar field 
interacts directly with matter particles through a confor- 
mal coupling of the form e^' Mpl . On the other hand, BD 
theory in the Einstein frame behaves like an interacting 
quintessence with j3 = — 1 / V6 + 4w 0] ■ More specifically, 
when the coupling constant is zero, the conformal cou- 
pling takes the form e -<f>/V6M pl i e ^ = _i/^/g. Thus, 
we expect that in the Einstein frame, the massive BD 
theory behaves as a Chameleon theory, however, is it 
true in the Jordan frame? It is worthy to stress that 
the physical results can be completely different in these 
frames and we should be very careful in interpreting the 
meaning of them. Although, it is not yet clear that which 
one of these frames are really physical but it is usually 
convenient to work in the Einstein frame which is math- 
ematically simpler and then convey the results to the 
Jordan frame. For more details see [2| chapter 2. 

As we mentioned before, the metric f(R) theory is 
dynamically equivalent to the massive BD theory. In 
the context of this theory, several authors claimed that 
an ill-defined behavior in the Newtonian limit exists 
(7 = 0.5) This conclusion can be inferred directly 
from analogy with BD theory (a; = 0). However, us- 
ing the Chameleon mechanism in the Einstein frame, 
Faulkner, T. et al [l5| have shown that the metric f{R) 
gravity can lead to 7 ~ 1. They have transformed their 
results to the Jordan frame. 

Here, using exactly the same method of fl6( | , we want 
to explore the Chameleon behavior and the thin-shell ef- 
fect in the Jordan frame of the massive BD theory with 
ui = 0. Assuming that the pressure is zero (p = 0), one 
can rewrite equation (|4|) as 

3D<j) + 2V - U-l)^- = R-k 2 p (19) 
dep 

Note that R = Now, if (f) ~ 1 and I X I < 1 then the 



above equation can be written as 

3Dcj) ~R- k 2 p (20) 

Thus, this theory will behaves like GR if 3 \0<j>\ < k 2 p. 
In the weak field limit, where one can use □ ~ V 2 , This 
condition can be expressed as follows 

A 2 | J?=fc2 ,V 2 p«p (21) 

where A^ = 1/rrio is the Compton wavelength of the field. 
This condition is the Compton condition. In order to in- 
terpret this condition, consider a source which its density 
changes slowly with radius say p ~ r £ (this is the case for 
the Earth and the Sun). Then the Compton condition 
can be written as 

A | fl=fe 2 p ^< P (22) 

Thus, if Ar is a length scale in which the density 
changes concretely (Ap ~ p) then we can infer that the 
density changes on the length scales that are much longer 
than the Compton wavelength. Physically it means that 
the scalar field interaction will suppress on scales larger 
than A^. 

As in [lq. in order to find the corresponding thin-shell 
condition [7[, we take the static spherically symmetric 
metric around a given source at the origin as follows 

ds 2 = -[1 - 2 A{r) + 2B{r)]dt 2 + [1 + 2 A{r)]{dr 2 +r 2 dQ) 

(23) 

We will assume that |A(r)| <C 1 and |-B(r)| -C 1 near 
the source and inside it. By taking into account these 
assumptions and the definition of the Ricci tensor, one 
can verify that 

V 2 (A + B)~-^R (24) 

V 2 B^- l -(R» + ^) (25) 

By substituting R® from 0-0 component of equation ([3]) 
into ([2~5"1) . we obtain 

(26) 

On the other hand, by using the conditions <fi ~ 1 and 
|^| < 1, we get 

V 2 B(r)^-\{R~k 2 p), 

I 1 (27) 

\7 2 A(r)c--(R-k 2 p)--k 2 p. 

And the field equation of the scalar field takes the form 

\7 2 0~±(R-k 2 p) = -\7 2 B(r) (28) 
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Thus, assuming that B(r) is finite at the origin, the so- 
lution of B(r) and the scalar field are related as 



B(r) = <fo - cb(r) = -A0(r) 



(29) 



where, as before. <j>o is the background value of the scalar 
filed. Following [16J , it is convenient to define an effective 
mass as 



l eff 



p{r') 



R(r') 



dv' 



(30) 



Thus, the solutions of A and B become 

2G [W)~np-W 2Gm eff 



B(r) 



A(r) 



\r — r | 
Gm Gm e f f 
r 3r 



(31) 



important to stress that although in the Jordan frame 
there is not any coupling between the matter and the 
scalar field but the scalar field couples directly to the 
Ricci curvature and consequently mass of the scalar field 
can depend to the curvature. Thus, we expect something 
similar to Chameleon effect in the Jordan frame. How- 
ever, this " Chameleon- like" behavior ensures us that the 
massive BD theories with vanishing coupling constant (or 
equivalently metric f(R) theory) can be compatible with 
the solar system tests. Also, in the cosmological con- 
siderations it can be different and distinguishable from 
ACDM. For example, phantom division can occur in it 

0- 

Now, we are in a position to answer the main ques- 
tion of our paper. Is it allowable to linearize the scalar 
field equation when the Chameleon-like behavior exists? 
From the thin-shell condition and the PN expansion of 
the scalar field we have 



where m = J p(r')dv' . Finally, the PPN parameter 7 is 
found to be 



.4 



7 : 



3m + \rn e ff 
A — B 3m + m e ff 



(32) 



It is clear form this equation that if m e // <C m then 
7 — > 1. An important result may be obtained from (|31[) 



l*(r)| <™j 



2G [PVW = 2 (33) 



Where $iv(r) is the Newtonian potential. Now, using 
equation (|2"9")) we get 



\A<f>\ < -§ N (r) 



(34) 



This condition sets an upper bound on the difference 
between the values of the scaler field from the interior to 
the exterior of the source. Note that if |A</>| <C |$^r(r) 
then m e f f <C m, so R — > k 2 p and we will recover the GR 
results. This condition is called this the thin-shell con- 
dition. It is easy to convert this condition to the corre- 
sponding Einstein frame version. Then it will be exactly 
the thin-shell condition obtained in with f3 = -1/V6- 
As has been discussed in [3, EH] , if this condition satis- 
fied then the exterior field is only generated by the thin- 
shell near the surface of the spherical source. However, 
it is worthy to note that we have named this behavior 
"Chameleon-like" here because it is different from what 
is known in the original Chameleon theories Q. In fact, 
in the Jordan frame there is no direct interaction between 
the matter and the scalar field and it looks like that we 
have strongly restricted the scalar field evolution such 
that it behaves like GR. Also, we will show that if the 
thin-shell condition satisfied then the scalar field needs 
to be heavy at the cosmological scales. But, we know 
that in the original Chameleon theory Q, the thin-shell 
condition do not makes any restriction on the interaction 
range of the scalar field far away from the source. It is 



(35) 



On the other hand, For having the Chameleon-like be- 
havior, ~ 1, hence by using (|35p we get 



< 1 



(36) 



Thus, as discussed before, the PN limit is applicable. 
In the following, we show that the second necessary con- 
dition for "Chameleon-like" behavior , i.e. |V'/V| ^> 1, 
forces the scalar field to be heavy at the cosmological 
scales. To do this We write the Ricci scalar as follows 



R(r) 



dV 



(37) 



Therefore, the background curvature scalar is Ro = Vq. 
Regarding that cf>o is the minimum of the effective poten- 
tial V e ff and also defining ?/(</>) and M 2 as 



v{<t>) 



M 2 = 



(38) 



we have 



M 2 
3M 2 



2Vr 



(3) 



(3) 

110 = ~14 



12 ism 2 



WVr 
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3M 4 2V r 



(-1) 



(4) 



48 120A'/ 2 120K 
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15M 4 



20V r 
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"0 



10M 2 F (3) 
Wo 



2V 
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Expanding ?i(0) in powers of A0 and using (|39|) . we get 
V' 2 



r]{4>) = 



M 2 



V O 



A0 



2Vp (3) A0 2 
W 2! 



(4) 



3M 4 A0 3 



(40) 



(5) 



10M 2 V r „ (3) 



A0 4 



O(A0 5 ) 



Note that we have expanded ??(0) to the higher orders 
of A0 because when the Chameleon-like behavior occurs 
the difference between |?/(0)| 3> 1 and r) ~ 2 is very 
large. First assume that M 2 <C 1, then one can rewrite 

rj((j)) as 



-(1 



K (3) Acj) 2 



V {4) A0 3 



K (5) A0 4 



v ' 



2! y ' 3! 
For \rj(4>)\ » 1 it is necessary that 



4! 



V 



(3) 



A0 2 
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V (4) A0 3 
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V, 
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A0 4 



» 1 



(41) 



(42) 



There are many cases for which this inequality is satisfied. 
For example, assume that these terms are of the same 
order, i.e. 



V, 



(4) 



(3) 
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(A0- 



> 1 



(43) 



Thus it is ,in principle, possible to satisfy |?y(0)| 1 
while ttiq/Hq 1. However, some notes are in order 
here. 

Does tuq/H 2 <C 1 guaranty that the scalar field is long 
range? Consider the field equation of the scalar field (Ql . 
We rewrite it as follows 



□ 0-m 2 (0-0 O ) + 



00 Vr 
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k 2 p 
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(44) 



If the second term is larger than the self-interacting 
terms, then the solution of for an isolated gravitat- 
ing system contains the Yukawa-like term e~ m ° r /r. So, 
the magnitude of mo is directly related to the range of 
the scalar field in the sense that if the scalar field is light 
(heavy) then it is longc (short) range. On the other hand, 
if the second term is negligible compared with the third 
term, even for small A0 , then the dynamics of the scalar 
field is dominated by the third term and there is not nec- 
essarily a Yukawa-like solution containing nio in it. In the 
other words, in this case, the lightness of the scalar field 
dose not necessarily mean the long rangeness of it. But 
for the Chameleon scalar fields Q, the lightness (heavy- 
ness) of the scalar field implies the long (short) rangeness 
of it. And of course, this is one of the main goals of that 
theory. 



The second problem for the Chamelcon-like behavior 
by the light scalar fields is that we encounter with non- 
smooth varying potentials, see (|4*3"|) . Thus, it is very 
unlikely to produce the Chameleon-like behavior by the 
light scalar fields (rn^jH 2 <1). 

Now, assume that M 2 ^> 1. In this case one can 
rewrite (|4U|) as 



, ± . 2 M 2 2 / V o (3) A0 2 
0o 2 0o V 2! 



y (4) A0 3 
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)- 
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It is clear that without any restrective conditions such 
as (|4"3"|) on the form of the potential and only with the 
assumption that itiq/Hq > lwe can produce |?7(0)| ^> 1. 
In this case one can write the Ricci scalar as follows 



Rq 



AH 2 



^( 2 )(r) + O(A0 2 



(46) 



Where Hq is the Hubble parameter. Experimentally 
we know that ^ = R ~^f" ^ 1 in the solar sys- 
tem (it is obvious in GR). For example on the Earth 
^ ~ OipAir/H^M 2 ^ ~ 10 27 and in the solar system 
we can estimate this ratio using the local Dark matter 
density which yields ^ - 0(p DM /H§M^) ~ 10 6 . Con- 
sequently, we can again infer from (|46j) that 



m| 
Hi 



> 1 



(47) 



which can be consistent with the condition (fTTj) . This 
condition shows that the Compton wavelength is very 
smaller than the Hubble length scale. This means that 
the scalar field is heavy in the cosmological scales. For 
example, near the Earth AR/R^ ~ 10 27 . By assuming 
that L rsj REarth and ip^ ~ 10~ 14 (for this assumption 
we have used (pJB]) ) and using equations (j3"5)l and (|4"6"|) it is 
easy to show that itiqL 2 ~ 10 3 . It is clear from (|4"6"|) that 
if the scalar field is very light in the cosmological scales 
then a gross violation of experiment will occur. We saw 
this result also in the PN limit where 7 was equal to 0.5 
if m\L 2 < 1. 

Thus, it seems like that the Chameleon-like behavior 
can olny occur for the heavy scalar fields (m^/H 2 3> 
1). And this is completely against the Chameleom scalar 
fields 0. As we mentioned before, the Chameleon scalar 
fields arc light and free at the cosmological scales. 

Also since 0o is the minimum of the effective potential, 
so we can write 



V 



"0 



V ' 



k 2 p 



(48) 



Where po ~ 10 29 h 2 is the current matter density of 
the universe. By multiplying this equation with L\ , Ll 
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is a length scale of the same order of the solar system i.e 
Ll ~ 10 15 cm, we can rewrite (|4"5)) as 

%l « 1 (49) 

which is equivalent to the condition (|12p which came from 
the PN limit considerations. 

So, we see that the trace of Chameleon-like behaviour 
is clear in the PN limit as we expected. As a final result, 
the PN limit of massive BD theory (with o> = 0) is appli- 
cable for the gravitational systems such as solar system 
where the Chamclcon-like behaviour can occur. And one 
can linearise the field equations to the PN order and the 
result of this expansion and the Chameleon behaviour is 
similar. As we mentioned before the Chameleon-like be- 
haviour which exists in the Jordan frame is different from 
the original Chameleon theory 0. On the other hand, 
we showed that the required conditions for this mecha- 
nism and those coming from PN limit considerations for 
having a viable theory in the solar system are consistent 
with each other. 



IV. CONCLUSION 

In this paper wc have shown that although the field 
equation of the scalar field in the massive BD theory 
allows the existence of the odd-order terms in the PN 
expansion of the scalar filed but they are all zero, this 
means that there is no non-linear terms in the field equa- 
tions. This is an interesting result because unlike the 
claim existed in the literature one can linearise field equa- 
tion to the PN order and also see the Chameleon-like 
behaviour. It is worth to mention that our conclusion 
dose not mean that the non-linearity in the expansion of 
the Ricci scalar, coming from the largeness of the local 
Ricci scalar relative to the cosmological background Ricci 
scalar, is not important. This means that smallness of 
y>( 2 ) relative to </>o dose not necessarily imply AR/Rq <C 1 
and one should note that the condition AR/Rq ^ f 
has been saved in the obtaining of the PN approxima- 
tion. For the massive BD theory for which the coupling 
constant is zero, we showed that the Chameleon-like be- 
haviour can exist if (f> ~ 1, |V/V'| <C 1 and if the thin- 
shell condition \ <€. ®N(r) satisfied and more impor- 



tantly this mechanism is different from the Chameleon 
mechanism. In fact, if the scalar field is forced to be 
heavy at the solar system through the Chameleon-like 
behavior then it has to be also heavy at the cosmological 
scales. And this is compeletly against the nature of the 
Chameleon scalar fields. Also, we showed that the results 
of PN limit and Chameleon-like behaviour are similar and 
when the Chamclcon-like behaviour exists, the perturba- 
tion of 4> due to the local system is very smaller than the 
cosmological background value. Thus, we can linearise 
the field equations to the PN order. 

We finish the conclusion by two example which con- 
firm the consistency of the PN consideration and the 
Chameleon-like behaviour. Consider the model f{R) ~ 
R 1+s . This model can be consistent with the solar system 

observations if 5 1.23 T 2.05 x 1CT 17 p} (Albeit we 

have taken into account the term ( 1 — 1 / vT — e 2 ) ~ 1 which 
is absent in the equation (83) of [l7j|) thus Chameleon 
mechanism occurs for this model. On the other hand 
from, (|46[) and the amount of AR/Rg in the solar sys- 
tem, we can write 

^| > 10 12 (50) 

By using this condition for the above model we reach to 
5 < 10~ 12 . 5 ~ 10~ 12 is sufficient for being consistent 
with the light deflection experiments, but for consistency 
with other observations such as the perihelion preces- 
sion of Mercury, smaller amount of 5 is needed. So, the 
condition (|50|) leads to a right bound on S in agreement 
with the result of [l^l • For another example consider the 
model f(R) ~ R - Ai? c (-^) P [H|. The bound on p com- 
ing from the Chameleon effect is p < 10~ 10 ^. On the 
other hand, using equation (|50[) one can easily verify that 
p < 10 -12 which is consistent with the previous bound. 
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